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1. Introduction 

In 1982 Thorne and MacDonald0 discussed the integral formulation of Max- 
well's equations in curved spacetime and used it in the following years to develop 
a number of applications to black hole physics, culminating in the "membrane 
paradigm" for treating the black hole horizona and physics around black holes. 
Their approach to the integral version of Maxwell's £auations is based on the 
Ellis 1 + 3 congruence description of those equationsBa but relies on the black 
hole spacetime splitting by the hypersurface orthogonal locally nonrotating ob- 
servers, also called the zero angular momentum observers (ZAMOs). Here this 
approach is generalized to an arbitrary observer family in any curved spacetime 
using the framework of gravitoelectromagnetism,cl simplifying some of their cal- 
culations by using an elegant differential-form transport theorem to calculate 
the time derivative of surface integrals. With this machinery, the meaning of 
the lines of force for observer-dependent electric and magnetic jields is clarified. 

A partial step in this direction was taken by Van BladelJH who developed 
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an extensive body of relativistic applications of Maxwell's equations in moving 
reference systems in Minkowski spacetime with an eye towards realistic engi- 
neering problems, with some limited steps towards extending them to curved 
spacetime, Jieavily influenced by the 1 + 3 discussion of Landau and Lifshitm 
and M0llerLl on Maxwell's equations in stationary spacetimes. The classic ex- 
ample of Maxwell's equations in a rotating coordinate system in Minkowski 
spacetimeB as well as static-observer measured electromagnetic fields in general 
stationary axisymmetric spacetimes both lead to integral Maxwell equations 
involving a pair of observer families. One is the general observer associated 
with the definitions of the electric and magnetic fields through the Lorentz 
force law on test particles and hence lines of force of the electric and magnetic 
fields, and the other is the hypersurface-forming observer associated with the 
hypersurface-submanifolds over which the integration is performed and used to 
define the flux integrals. When these observers are not the same, the usual 
flat-spacetime inertial-coordinate relationships between lines of force, flux and 
sources no longer hold. 

In this article differential form language coupled with a modern version of 
the various integral transport theorems is used to obtain the set of generalized 
integral Maxwell equations. Some observations about lines of force are then 
made in concluding remarks. An appendix shows the relationship of the modern 
statement of the transport theorems to those from vector analysis. 

2. Maxwell's equations in 1 -f 3 form 

The Maxwell equations can be expressed covariantly in many ways. In dif- 
ferential form language one has 

dF ^0 , d*F = -47r (or 5F = -inJ^ ) , (1) 

where F is the Faraday electromagnetic 2-form and J is the current vector field, 
obeying the conservation law 

SJ^ = *d*J^ = . (2) 

Here 5 = *d*is the divergence operator, * is the duality operation on differential 
forms (or their associated index-raised contravariant tensors), and jj and b are the 
index raising and lowering operations taking a tensor to its fully contravariant 
or fully covariant form. 

The splitting of the electromagnetic 2-form F by any observer family (with 
unit 4- velocity vector field u: u ■ u = —1, metric signature — h ++) gives the 
associated electric and magnetic vector fields E{u) and B{u) as measured by 
those observers through the Lorentz force law on a test charge, and the relative 
charge p{u) and current density vector field J(it). This splitting procedure is 
a standard orthogonal decomposition of tensor fields associated with the or- 
thogonal direct sum of each tangent space into the directions along u and the 
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complementary directions in the local rest space of u. This "relative observer 
decomposition" of F and its dual 2-form *F is 

F = [uA E{u) + , 
*F = [-U A B{u) + *^^^E{u)f , 

while J has the representation 

J = p{u)u + J{u) . (3) 

Here is the spatial duality operation over the local rest space of the 
observer family. For example, in an observer- adapted orthonormal frame {ca} 
(a = 0, 1, 2, 3) adapted to u = eo, these definitions are equivalent to 

Eiur = F% , B{ur = - *F% , piu) = J" , j{ur = piur^j^ , (4) 

where 

P{u)'-(i = 5''p+u^U0 (5) 

is the spatial projection tensor. 

If f} is the unit oriented volume 4-form (770123 = 1 in an oriented orthonormal 

frame) and ri{u) = tt—lr; is its left contraction with u (the spatial volume 3-form), 
equivalently ri{u)aj3'y = u^VSaff-y, then the spatial cross product of two spatial 
vector fields is defined by the successive right contractions 

^x„r = {[r/(u)Lr]i_x}» = *(")[XAy] , (6) 

or in components 

[X XuYr = 7^{urp^X^Y'' , (7) 

where ?7(m)i23 = 1 in an oriented, spatially oriented observer- adapted orthonor- 
mal frame. The spatial duality operation is then 

Biur = Iviur^.F'^^ . (8) 
Similarly the spatial dot product of two spatial vector fields is 

X{u) •„ Y{u) = P{u)^pX{urY{u)f . (9) 

If U is the 4- velocity of any test particle with charge q and nonzero rest mass 
m, it has the orthogonal decomposition 

U = ^{U,u)[u + u{U,u)] , (10) 

defining its relative velocity and associated gamma factor with respect to u. The 
absolute derivative of U with respect to a proper time parametrization along 
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its world line is its 4-acceleration a(U) = DU/drij. The Lorentz force law then 
takes the form 

ma{U) = qj{U,u)[E{u) + iy{U,u) x„B(7i)] . (11) 

Using the notation V{u) = P(u)V for the spatial projection of the covariant 
derivative, namely 

[y{u)xYr = X'P{uysP{urpV^YP , (12) 

the orthogonal splitting of the covariant derivative defines a time derivative 
V(fw) (u) = P(m) V„ (the spatial Fermi- Walker derivative) and a spatial covariant 
derivative vector operator V{u) such that V(u) = V(u)x- One then has the 
generalizations of the usual spatial vector analysis operators grad.^j/ — V(u)/, 
curl„X = V(u) x„ X, and div^X — V(u) •„ X. The spatially projected Lie 
derivative along u is also needed and is given the analogous symbol £„(u) = 
P{u)£u- These differential operators have been studied by Bini, Carini and 
JantzenB 

The relative observer formulation of Maxwell's equations is well known. Pro- 
jection of the differential form equations (||) along and orthogonal to u gives the 
spatial scalar (divergence) and spatial vector (curl) equations. Using the nota- 
tion and conventions of Bini, Carini and JantzenEl these are 

div.aB{u) + H{u) -u E{u) = , 

curluE{u) - g{u) x„ E(u) + [£{u)u + e{u)]B{u) = , 
div„£;(M) - H(u) -u B{u) = Anp{u) , 

curl„S(u) - g{u) x„ B{u) - [£(u)„ + Q{u)]E{u) = AitJ{u) , (13) 

where the gravitoelectric vector field g{u) = —a{u) = — V„m and the grav- 
itomagnetic vector field H{u) = 2[*Wuj{u)^Y of the observer u (sign-reversed 
acceleration and twice the vorticity vector field) are defined by the exterior 
derivative of u 

du^ = [uA g{u) + *^^m{u)f . (14) 

The expansion scalar Q{u) = —6u^ = Ti 9{u) appears in an additional term in 
the covariant derivative of u as the trace of the (mixed) expansion tensor 0{u), 
of which the shear tensor (t{u) — 6{u) ~ ^Q{u)5 is its tracefree part 

Vu = -a{u) ®v° + 0{u) - uj{u) , 

u";0 = -a{u)''uf3 + e{ufp - lo{uTp . (15) 

The sign of the (mixed) vorticity tensor lo{u) depends on the index, signature 
and orientation conventions, and is chosen here so that its spatial dual is the 
commonly accepted vorticity vector; in an oriented frame one has 

u:{ur = \usTi'''^''u:{u)p^ = \sif'^P''u^,p = [icurl„7.]" . (16) 
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J his representation of Maxwell's equations differs from the Ellis representa- 
O only in the use of the spatially projected Lie derivative rather than the 
spatially projected covariant derivative along u (spatial Fermi- Walker deriva- 
tive) . These two derivative operators are related by the following identity for a 
spatial vector field X (orthogonal to u) 

[£iu)u + e{u)]X = [V(fw) (u) + {-a{u) + u;{u)}l.]X . (17) 

Black hole spacetimes have two privileged observer families, the static ob- 
servers with 4- velocity m satisfying 0(m) = = cr(rn), and the zero-angular- 
momentum observers (ZAMO's) with 4-velocity n satisfying 0(n) = = tj(n). 
For both u = 71,111, there exist acceleration potentials (lapse bmctions) 

g{u) = -grad,„lnL(u) , L{n) = N , L{m) = M , (18) 

which allow the curl and acceleration terms in Maxwell's equations to be rewrit- 
ten as a single term 

L{uy^cur\u[L{u)X{u)] = curl„X(w) - g{u) Xu X{u) . (19) 

In a comoving coordinate system {i, a;*} (i — 1,2,3) adapted to u, then the 
observer four- velocities and Lie derivatives simplify to 

n = N-'[d/dt - N] , £{n)„ ^ N'^ £{n)^^g^_^ , 

m = M-'[d/dt] , £{m)m = M'^£im)Q/9t , (20) 

where is the shift vector field. 



3. Spacetime integrals in 1 + 3 form 

The integral form of Maxwell's equations in classical physics relates the 
fiux of electric/magnetic field through a closed surface in space to the enclosed 
electric/magnetic charge and relates the time derivative of the flux through a 
loop in space to the line integral of the magnetic/electric held around the loop 
and the current passing through the loop. These integral equations are obtained 
from the four differential Maxwell equations by applying Gauss's divergence 
theorem to the two scalar equations and Stokes' theorem to the two vector 
equations. 

Developing the analogs of these equations in a curved spacetime requires an 
observer family to split the electromagnetic field and a slicing of the spacetime 
by a family of timelike hypersurfaces, which may be taken as the hypersurfaces 
associated with constant values of some time function t. The latter structure 
is needed to describe the curves, surfaces and spatial regions over which the 
integration is performed at some moment of "time." If the observer family is 
taken to be the field n = —Ndt'^ of imit normals to the slicing (4-velocity field 
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of the "hypersurface observers"), then much simpHfication occurs, but if the 
desired observer family is not hypersurface orthogonal, as is the case even for 
rotating reference frames in Minkowski spacetime, one is forced to consider the 
more general case. Over the years even this simple case has been the subject 
of many confusing discussions in the literature regarding Maxwell's equationsO 
and even relativistically "correct" rotating frame transformations. 

In the integral formulation of Maxwell's equations, Stokes' theorem is applied 
to the 3-form exterior derivatives of F and *F using in the spacetime differential 
form version of Maxwell's equations (|l|). Two distinct types of hypersurfaces 
must be used in this application in order to yield the spatial scalar and spatial 
vector information with respect to a time function and a test observer family. 
Let V{t) denote a family of spatial regions parametrized by the time t, filling 
a world tube in spacetime, with boundaries dV{t) which are closed 2-surfaces 
at time t with the outward orientation. This first type of hypersurface leads to 
the integrated form of the spatial divergence equations projected along the unit 
normal n, even if expressed in terms of a different test observer family. 

The second type of hypersurface is the world sheet of a family of 2-surfaces 
S{t) between two hypersurfaces t and t + e. The boundary of this world sheet 
consists of S{t + e) and S{t) with opposite inner orientations. Taking the limit 
of Stokes' theorem here as e — > leads to the time derivative integral equations, 
with the integral on the 2-surface composed of the family of boundaries dS{t) 
between t and t + e reducing to a line integral on dS{t) alone. These will depend 
on the relative velocity of the 2-surface and the unit normal n. An alternative 
to this limiting operation is to use a transport theorem. 

To discuss these questions uniformly, let P'-^^ (t) denote a family of p-surface 
integration domains within the family of time slices and parametrized by the 
time t, namely a family of spatial regions V{t) (p—i) or a family of surfaces S{t) 
(p=2) or a family of closed curves C{t) (p=l). Let dC{t) denote the correspond- 
ing family of boundary (p— l)-surfaces. Let U be an evolution vector field whose 
1-parameter family of diffeomorphisms drags the family 2?*^^^ (t) into itself, satis- 
fying £u t = U Adt = 1. Dragging along by U is then equivalent to translation 
in the coordinate t when completed to a system of comoving coordinates {t, a;*}. 
Assuming that U is timelike (or at least nonnuU), the corresponding 4- velocity 
(or unit tangent) is 

U ^-f{U,ri)[n + v{U,n)]=-f{U,u)[u + v{U,u)] , (21) 

so that 

U ^N-i{U,n)-^ij , (22) 

where N = —n ■ U is the slicing lapse function. Such evolution vector fields U 
can differ by vectors which are tangent to the V'^p'^ {t), corresponding to an extra 
motion of the region or surface or curve into itself. If Ut denotes the 1-parameter 
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family of diffeomorphisms associated with U (its "flow"), then 

D(f)(t) = C/t_t„pW(to) . 



(23) 



A similar equation holds for the corresponding boundaries. 

For a spacetime p-form O and p-surface integration domain V^^^ {t) within 
the hypersurface of time t, the transport theorem is then 



d_ 
di 



f 



where 



£u^ = UAdVL + d[UAn] 



(24) 



(25) 



is a useful formula for the Lie derivative of any p-form, and which allows the 
second term to be integrated to the boundary using Stokes' theorem 



av(p){t) 



C 



(26) 



where ^ is a (p — l)-form. 

The differential form integrals require no metric, but they can be converted 
to more familiar but equivalent metric notation for volume, surface and line 
integrals, which is related to the p-foim integral notation for a function / and 
a spatial vector field X{u) by 



/ fv{u) = f fdV{u), 

Jv{t) Jv(t) 

[ *(«)X(u)^ = / X{u)-udS{u) 
Js{t) Js{t) 

[ X{uf = [ X{u) -u di{u) , 

Jc(t) Jc(t) 



(27) 



IC{t) JC(t) 

For example. Stokes' theorem for p = 2, 3 then becomes 



/ curl„X(w)-„rf5(w) = / X{u)-udt{u), 

Jsit) Jas(t) 

[ divu X{u)dV{u) = [ X{u)-udS{u). 

Jv(t) JdV(t) 



/V(t) JdV{t) 

The transport theorem is the result of the calculation 



(28) 



dt 



dt 
d_ 

~ dt 



t=to JUt-toV(p>{t„) 



[ UiXn= [ £un 



(29) 



Integral formulation of Maxwell's equations 8 



transforming the integral and using the definition of the Lie derivative. This 
single transport theorem corresponds to the line, surface and volume integral 
formulas found in some older textbooks on electromagnetism and collected to- 
gether and expressed in conventional notation in appendix D of van BladelQ and 
discussed in the appendix. 

4. Integral 1 + 3 form of Maxwell's equations 

The "spatial scalar" Maxwell equations in integral form come from integrat- 
ing Maxwell's equations (0) over a 3-surface region V{t) of a time coordinate 
hypersurface. Integrating first dF = leads to 

0=/ dF = F^ [uA E{u) + *i"W{u)f . 

Jv(t) JdV(t) JdV(t) 

However, restricting the 1-form to a constant time coordinate hypersurface 
implies 

= -Ndt\t=to = n^\t=ta = l[n,u)[v:' + }J'{n,u)\\t=ta ^ "1t=to = -^^^ {n,u)\t=to 

(30) 

so the right hand side of the integral equation becomes 

dV{t) 

We then have 

• Gauss's law for the magnetic field 



/av(t)[-^(") - ^("'") ^-^H = • (31) 

The final expression uses the more conventional metric notation for volume, 
surface and line integrals. 

Analogously, the integral Maxwell equation dual to this one comes from 
integrating the dual Maxwell equation. The right hand side of this integral is 

-An [ *f = -47r ( *[p{u)u + J{u)]^ 

Jv{t) JV{t) 

= An *(")[p(u) - i/(n, u) •„ J{u)] 
Jv{t) 



= An I [p{u) — v{n,u) •„ J(u)] dV{u) , 
Jv{t) 

so that one has instead (the spacetime duality operation on F is equivalent to 

[E{u),B{u)] ^ [-B{u),E{u)]) 

• Gauss's law for the electric field (Coulomb's law) 
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[E{u) + iy{n, u) x„ B{u)] •„ dS{u) 
dv(t) 



= 47r / [p{u) - iy{n, u) •„ J{u)] dV{u) . (32) 
Jvit) 

The extra velocity terms in both Gauss law equations correct for the motion 
of the observers with respect to the time hypersurface in which the integration 
takes place. 

The "spatial vector" Maxwell equations in integral form come from applying 
the transport theorem ( p^ ) with fl = F and ^ *F respectively on a 2-surface 
S{t). In the first case, the left hand side of this equation is 



"I Js(t) "t Jgu) 



5(t) JS{t) 

lb 



d 

di 
d 
di 



x„ E{u) + B{u)Y 

S(t) 

[B{u) ~ v{n,u) Xu E{u)] dS{u) , 

S(t) 



while the right hand side is 

/ £uF^ I [U -IdF + d{U -IF)] ^ I U-IF, 

Js{t) Jsit) JdS{t) 

where Maxwell's first equation and Stokes' theorem have been used. Using 
and (El]), this 1-form evaluates to 



= N-f{U,n)-'^U_iF = ~Nj{U,n)-^E{U) 

= -N-/{U,n)-'^-f{U,u)P{u,(j)-^[E{u) + iy{U,u) x.^B{u)] 

= -N'y{U,n)-^-f{U,u)[E{u) + v{U ,u) x„ B{u) + u[v{U ,u) ■ E{u)]]^ 

so that 

U^F\t^to = ~N-i{U,n)~^-i{lJ,u)[E{u) + v{U,u) x^Biu) 
-:^{n,uMU,u)-E{u)]f \t^to , 

where the "inverse" of the mixed projector P{u, U) — P{u)P{U) and the electric 
field transformation law have been given by Jantzen, Carini and Binicl (equations 
(4.7) and (4.14)) and equation ( ^0| ) has been used. The final form is then 



• Faraday's law 

[B{u) — i^(n, u) Xu E{u)] -u dS{u) 



d_ 

' dt 



S(t) 

Nj{U, nr^-f{U, u)[E{u) + v{U, u) Xu B{u) 

dS{t) 

-v{n, u)[v{U, u) ■ E{u)]] -u dl{u) , (33) 
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where the right hand side merely represents the hnc integral of the component 
of the Lorentz force along the direction of the curve, namely 

N[E{n) + v{U, n) x„ B{n)] •„ dl{n) , (34) 

corrected for the time t rather than the proper time associated with n and then 
re-expressed in terms of the generic observer u. In fact one can write this right 
hand side in the alternative form 

L{u)[E{u)+U x„ {B{u) - v{n,u) XuE(u)}] •„ d^(u) , (35) 

ds{t) 

where the "lapse-like function" L{u) ~ N^{u, n)~^ reparametrizes the u proper 
time derivative to the t time derivative in this equation and 

7([/,n) l{U,n) 

is the 4- velocity of the curve dS{t) reparametrized by the u proper time and 
projected into the local rest space of u, which would have spatial coordinate 
components W = dx^ / dTu in comoving coordinates adapted to u (denoted by 
V in Thorne and MacdonaldB). The extra term containing it corrects for the 
motion of the curve with respect to the observer u. 

The dual equation is obtained in a similar way, with the additional term 
resulting from Maxwell's second equation 

UAd*F = -47r / *f 
Sit) Js{t) 

= -Att [ Nj{U,ny^j{U,u)[-p{u)i^{U,u) + J{u) 
Js{t) 

+1/(71, -u) Xu [iy{U,u) Xu J{u)])] -u dS{u) , 

where the relations 

C/_l = C/_l - (C^)J(C/)^ (37) 

and 

ri{U) = 7(J7, u)[7]{u) -uA *(">;/([/, u)] 
have been used. This term can also be rewritten as above 

L{u) [J{u) - U{p{u) - iy{n, u) ■„ J{u)}] •„ dS{u) , (38) 

Sit) 

where again the U term corrects for the motion of the surface as above. 
The final form in the original notation is then 
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• Ampere's law 

d f 

-— [E{u) + v{n, u) x„ B{u)] dS{u) 
Js(t) 

Nj{U, ny^j{U, u) [B{u) - u{U, u) x„ E[u) 

dS{t) 

—v{n, u)[iy{U, u) ■ B{u)]] •„ dl{u) 
+47r / Nj{U,n)-^j{U,u)[-~p{u)iy{U,u) + J{u) 

u) x„ u) Xu J{u)]] -u dS{u) . (39) 

The integral form of charge conservation is obtained by applying the p — 3 
version of the theorem (29), integrating = on a hypersurface region 
V{t), using d*J^ = to eliminate the first term in the Lie derivative and then 
applying Stokes' theorem to the second term dU _l and finally using ( p7| ) 
to re-express U _l . Alternatively one can use Coulomb's law for the electric 
field to replace the flux integral by the charge integral in the preceding equation, 
where S{t) = dV{t) so that dS{t) is the empty set. The result is 

• Charge conservation 



dt 



^ / [p{u) - vin, u) •„ J{u)]dV{u) 

V(t) 

Nj{U, n)-^j{Lf, u)[-~p{u)iy{U, u) + J{u) 

dvit) 

+v{n, u) Xu [v{ij, u) x„ J{u)]] -u dS{u) . (40) 

It is important to note that each of the integrals which appear in the integral 
form of Maxwell's equations are merely expressed in terms of the observer with 
4- velocity u, but the values of these integrals are independent of this observer, 
and reduce to the values assumed when the observer family is the family of 
normal observers (m = n so v{n^u) — 0). These values are the Thorne and 
Macdonald equations^ 

B{n) •„ dS{n) = , 

dV(t) 

E{n) •„ dS{n) = 47r / p{n) dV{n) , 
dV{t) Jv{t) 

N[E{n) + iy{U, n) x„ B{n)] •„ di{n) = / B{n) •„ dS{n) , 

ds(t) "I Js{t) 

N[B{n) ~ iy{U, n) x„ E{n)] •„ diin) 

dS{t) 

4 / E{n) •„ dSin) + 47r / N[J{n) - , n)p{n)] •„ dS{n) , (41) 
dt Js{t) Js{t) 
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and if the surface S{t) is also at rest with respect to these observers {I'iU, n) — 
0), the latter two simplify further. If the surface S{t) is at rest with respect to 
the observers u, then U = u and ^{U , u) = 0, leading to a different simplication. 

5. Lines of force 

In Minkowski spacetime, it is customary to picture electric and magnetic 
fields at some moment of time by drawing their corresponding lines of force on 
space or some cross-sectional plane surface within space. The density of field 
lines (in charge-free regions for the electric field) can be taken to directly reflect 
the magnitude of the field strength since the integrated Maxwell divergence 
equations show that the flux is conserved, i.e., in any flux tube whose lateral 
sides (by definition) consist of field lines, the number of field lines through the 
ends of the flux tube is the same at both ends. Field lines are not created 
or destroyed. Thus the number of field lines passing through a fixed area can 
be taken as proportional to the field intensity. For a closed surface, the total 
electric flux directly determines the enclosed charge, while the total magnetic 
flux is zero. 

The same is true for the hypersurface observers in a general spacetime, since 
the same integral form of these Maxwell equations holds. The flux of electric 
and magnetic field lines is conserved and the integral of the charge density over 
a volume determines the enclosed charge, which is proportional to the electric 
flux out of its boundary, while the total magnetic flux out of a volume is always 
zero. The electric lines of force of a point charge at rest near a black hole and 
the magnetic lines of force around a black hole immersed in an asymptotically 
uniform magnetic field, as seen bv the ZAMO hypersurface observers, were first 
introduced by Hanni and RuffiniE3 and used by Rufhni and collaboratorsEj'oli3 
to stiidy electrodynamic properties of black holes (see also later work by Thorne 
et afl). 

However, for general observers, the separate fluxes of the electric and mag- 
netic fields are no longer conserved and one loses the nice connection between 
field intensity and the density of flux lines. Similarly the classic relation between 
electric flux and charge is lost, and the total magnetic flux out of a closed sur- 
face is in general not zero. Even the integral of the charge density as observed 
by u over a volume does not directly determine the enclosed charge. Instead 
one is forced to work with the combinations of the fields which are equivalent 
to transforming back to the hypersurface observer quantities, which is why the 
integral form of Maxwell's equations is so complicated for a general observer. 
However, van BladelQ has shown that in spite of this complication, it is still pos- 
sible to use them in describing certain kinds of problems in accelerated frames 
of reference. 

Acknowledgements 



Integral formulation of Maxwell's equations 13 



The authors thank Remo Ruffini for useful discussion and support. 
References 

[I] K. S. Thorne and D. Macdonald, Mon. Not. R. Astr. Soc. 198, 339 (1982) [plus 
microfiche MN 198/1]. D. Macdonald and K. S. Thorne, Mon. Not. R. Astr. Soc. 
198, 345 (1982). 

[2] K. S. Thorne, R. H. Price and D. A. Macdonald, Black Holes: The Membrane 

Paradigm (Yale University Press: New Haven, 1986). 
[3] G. F. R. Ellis in Cargese Lectures in Physics, Vol. 6 ed. E. Schatzman (Gordon 

and Breach: New York, 1973). 
[4] G. F. R. Ellis and H. van Elst in Theoretical and Obser vational Cosrr ioloqy ed. 

Marc Lachieze-Rey, (Kluwer: Dordrecht, 1999) pp. 1-116 fer-qc/981204^ ]. 
[5] R. T. Jantzen, P. Carini and D. Bini, Ann. Phys. (N.Y.) 215, 1 (1992). 
[6] J.F. Corum, J. Math. Phys. 18, 770 (1977); 21, 2360 (1980). 
[7] J. Van Bladel, Relativity and Engineering (Springer Verlag: New York, 1984). 
[8] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Permagon 

Press: New York, 1975). 
[9] C. M0ller, The Theory of Relativity, Second Edition (Oxford University Press: 

Oxford, 1972). 

[10] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (Freeman: San Fran- 
cisco, 1973). 

[II] D. Bini, P. Carini, R. T. Jantzen Class, and Quantum Crav. 12, 2549 (1995). 
[12] R. Hanni and R. Ruffini, Phys. Rev. D 8, 3259 (1973). 

[13] R. Hanni and R. Ruffini, Lett. Nuovo Cim. 15,) 189 (1976). 

[14] R. Ruffini, "Physics Outside the Horizon of a Black Hole" in Physics and As- 
trophysics of Neutron Stars and Black Holes, Eds. R. Ruffini and R. Giacconi, 
(North Holland: Amsterdam, 1977). 

[15] T. Damour, R.S. Hanni, R. Ruffini and J.R. Wilson, Phys. Rev. 17, 1518 (1978). 



Integral formulation of Maxwell's equations 14 



Appendix A. Time derivative integral formulas: transport theorems 

The elegant differential forms statement of Stokes' theorem simultaneously 
captures the vector analysis integral theorems for line integrals, surface inte- 
grals and volume integrals and their boundaries on a 3-dimensional Riemannian 
manifold. In a similar way the modern statement of the transport theorem con- 
tains the individual results in vector analysis language. First the flat spacetime 
versions of these are recovered, and then the hypersurface observer versions are 
obtained. 



5.1. Van Bladel integral formulas 



Some older textbooks on electrodynamics contain these transport theorems, 
restated by Van Bladeli. Starting from the the general form (|2^) of the transport 
theorem, specialize to flat spacetime and Minkowski coordinates {i, x*}, with 
n = dt usual unit time direction vector field, X{n) — X{nydi a spatial vector 
field, / a spacetime scalar, and the notation M(n) = dx^di^ dS{n) = 5'^^ tjkidx^ A 
dx^di, dV{ji) = dx^ A dx'^ A dx'^ . One finds 



- / X{n) •„ dl{n) 

JC{t) 



d 
di 



d 
di 



C{t) 

r 

S{t) 
Sit) 

f 

V(t) 
Vit) 



dX{n) 



dt 

X{n) •„ dS{n 



+ grad [X{n) ■ iy{n)] + [cm\X{n)] x i/(n) 



(A.l) 

diin) , 

(A.2) 



dXjn) 
dt 

fin) dV{n) 
df{n) 



+ i'{n)div X{n) + curl [X{n) x i^(ri)] 



dS{n) 



(A.3) 



dt 



where 



+ /(n)divi/(n) -f i>{n) ■ grad/(r7,) 



- + Kn)-grad/= — 



dV{n) 



(A.4) 



is the total derivative along the curve parametrized by coordinate time t and 
having the tangent vector U = dt + v{nydi. 

For example, consider deriving the second of these. Write the curl term in 
equation (A.2) as 

curl[A(n) x + i^(n)°[div A(n)] = £^(„) A(n)'' + A(n)"[div . (A.5) 

Taking into account the relation £Q^X{nY — dtX{n)°- and by adding this term 
to both sides of equation ( A.5 ) one finds 

dtX{nY + curl[A(n) x v{n)Y + v{nY[dw X{n)] 
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= £uX{n)'' + X{n)''[divi^{n)] . (A.6) 

Finally, contraction by dS{n)a — \'ri{n) abcdx^ A dx"^ gives equation ( [A.2| ) 

[dtX{nY + curl[X(n) x v{n)Y + v{nY[diY X{n)]] dS{n)a = 
= [£uX{nY + X{nY[A\Yv{n)]] ^f]{n)abcdx'' A dx" 

= *(")[£[/[77(n)l_X(n)]]b4(ix'' Adx^ , (A.7) 

for which 17(2) ^ 

Analogously equation (A^) corresponds to = f{n)rj{n). In fact one has 

£u[I{n)] = + v{n) ■ Wf{n) , (A.8) 

so that 

*^"^£u[fin)vin)] = + i^in) ■ V/H + (div^.(n)) /(n)] , (A.9) 

where the relation *(")£[/ [77(71)] = divj/(n) has been used. 



The first relation (A.l) corresponds to fi^^^ = X{n)'' . This is easily obtained 
starting from the vector identity 



[[curlX(72)] X iy{n)]^ = -t/(7i)^V„X(n), + iy{nYVcX{n)a 
and the Lie derivative relation 

pinrV,X{n)a = £,(n)X{n)a - X{n)i,Vav{nf , 

so that 

dtXa + [{cur\X{n)) x i.(7i)], = £uX{n)a - Va(X(n) • 1^(71)) 



(A.IO) 



(A.ll) 



(A.12) 



5.2. Thome and MacDonald 



The corresponding generalizations to any hypersurface observers in a general 
spacetime, where n is the unit normal to a slicing by spacelike hypersurfaces, 
are the Thorne-MacDonald equations (2.27), (2.26), and (2.25) respectively § 



d 
di 



d 



C(t) 



X{n) •„ di{n) 



(A.13) 



N 



C{t) 



V{n)nX{n) + e\_X{n) + [\J{n) x„ Xin)] x„ v{n) •„ dl{n) , 



dt 



- / X{n) •„ dSin) 



5(t) 



(A.14) 



N 



/5(t) 



y{n)nX{n) + [Sin) - dl.]X{n) + iy{n)[V{n) •„ X{n)] 
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+ V(n) x„ [NX{n) x i/(n)] •„ dS{n) 



d_ 
di 



f{n) •„ dVin) 



(A.15) 



V(t) 



v(0 



N[V{n)n + e(n)]/(n) + V(n) ■ [Nf{n)u{n)] dV{n) . 



Of course these equations reduce to Van Bladel's when N = 1, = = 
G(n). The precises Thorne-MacDonald form of the last two of these equations 
is obtained by using Stokes' theorem on the final curl and divergence terms to 
convert them respectively to a line integral and a surface integral. 



